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ABSTRACT 



The dynamic response of elastically-mounted cylinders in 
a sinusoidally oscillating two-dimensional flow has been 
investigated both theoretically and experimentally. 

The experiments were carried out in a large U-shaped water 
tunnel, with smooth and sand-roughened circular cylinders. 

The results have been expressed in terms of lift coefficients, 
relative amplitudes of oscillation, Keulegan-Carpenter number, 
Reynolds number, Strouhal number, and a mass parameter. 

A theoretical analysis has been carried out by assuming 
the transverse force varies harmonically. Furthermore, the 
instantaneous values of the displacement have been predicted 
through the use of the Duhamel's integral. 

The results have shown that elastically-mounted cylinders 
undergo synchronous oscillations at a reduced velocity of 5.7 
when the vortex shedding frequency coincides with the natural 
frequency of the cylinder. Away from the region of synchroni- 
zation the cylinder is subjected to forced response by vor- 
tices shedding at a frequency smaller or larger than the 
natural frequency of the cylinder. 
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Sum of the force transducer outputs (F^+F 2 ) 

Strouhal number (St = f D/U ) 

v m 

Fluid flow period 
Cylinder oscillation period 
Vortex shedding period 
time 

Duhamel superposition impulse time 

Fluid flow velocity 

Maximum fluid flow velocity 

Reduced velocity (U /f D) 

2 nr n 

Maximum cylinder displacement 

Maximum normalized cylinder displacement (Y^/D) 

Initial cylinder displacement 

Initial cylinder velocity 

Cylinder displacement 

Cylinder velocity 

Cylinder acceleration 

Velocity 

MC 

Mass parameter (A = — f—) 
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I. 



INTRODUCTION 



The investigation of the nonlinear deterministic or sta- 
tistical response of an offshore structure or of its tubular 
members is of considerable interest and may be carried out 
through mathematical modeling, small-scale experiments, and 
prototype testing or through a combination of these tech- 
niques . 

Much of what is known about the wave forces and moments 
acting on bodies and complex structures of various shapes 
came from field data and experiments in wave channels. To 
extrapolate from these to the prediction of forces on new 
structures in natural waves is often very uncertain, due to 
the difficulty or impossibility of modeling all parameters 
that might have an effect. For example, there is considerable 
uncertainty as to how to account for the effect of currents 
which may be superposed on waves / — 1_/, for the elasticity 
and flow induced oscillation of the members of the structure, 
for the roughness and diametral increase of the connecting 
members, for the effects of ambient flow with shear, ete . , 
not to mention the two most commonly discussed dimensionless 
parameters, namely, the Reynolds number and the relative 
amplitude . 

Studies of Reynolds number, relative amplitude, and 
roughness effects on the in-line and transverse force coeffi- 
cients have recently been carried out by Sarpkaya / — 2, 3 /. 
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These studies were prompted by the need for systematic data 
and partly by the controversy concerning the influence of 
Reynolds number and surface roughness on the time-dependent 
loads on circular cylinders in harmonic flow. In all the 
experiments, the cylinders were mounted rigidly enough so that 
there was no possibility of either in-line and/or transverse 
oscillation. 

The structural members in a prototype are seldom, if ever, 
designed to preclude the occurrence of forced or synchronized 
oscillations under the action of waves and/or currents. Evi- 
dently, the oscillating flow can cause vibrations in the in- 
line direction. More importantly, however, the transverse 
force, which is associated with vortex shedding, may give rise 
to serious transverse oscillations. The much-mentioned fail- 
ure / 4 / of Texas Tower No. 4, an offshore platform off the 
coast of New Jersey, on 15 January 1961, is a classic example 
of vibration of a marine structure induced by oscillating 
flow and currents. It is evident from the foregoing that 
methods of prediction must be developed to prevent catastro- 
phic flow-induced vibrations of the tubular members of marine 
structures subjected to wave action. 
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II. REVIEW OF PREVIOUS INVESTIGATIONS 



The non-linear vibration problems of elastically-mounted 
structures in oscillating flows have not been sufficiently 
explored. While much progress has been made regarding the 
in-line and transverse forces acting on stationary structures 
/~~ 2, 3_/ and regarding the oscillations of elastic structures 
in steady flows (see, e.g.. Refs / 5-1 _J .) , there has been 
relatively little theoretical and experimental work on the 
complex dynamic response of elastic structures to oscillating 
flows. This has been partly due to the lack of reliable data 
for the force transverse-coefficients for both stationary and 
vibrating structures . 

Laird / 8_/ explored, in 1962, the effects of support 
flexibility by oscillating a vertical cylinder through still 
water. He found that the forces acting on a flexibly-supported 
oscillating cylinder can exceed 4.5 times the drag force on 
the cylinder rigidly-supported while moving at a uniform 
velocity equal to the maximum velocity during the oscillation 
and that a cylinder, flexible enough to have transverse oscil- 
lations with amplitudes more than half the diameter, while 
performing large amplitude oscillations in water, tends to 
oscillate transversely at the eddy-shedding frequency and to 
vibrate at twice the eddy frequency in the in-line direction. 
Laird's work pointed out the significance of the vibrations 
in harmonic flow and the need for a detailed investigation of 
the phenomenon . 
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Vaicaitis / — 9_/ investigated the response of deep-water 
piles due to cross-flow forces generated by wind-induced 
ocean waves. The resulting cross-flow forces were treated as 
random processes in the time-space domain and were assumed to 
be dependent on fluid velocities and vortex shedding processes. 
Out of necessity, Vaicaitis had to make a number of assumptions 
regarding the structural motion, the added mass coefficient, 
Strouhal number, etc. Vaicaitis was careful to note that some 
of his conclusions were based on a particular example and that 
they might not be true for piles with different dynamic char- 
acteristics . 

Selna and Cho / 10_/ investigated the in-line resonant 
response of a tall structure by assuming constant drag and 
inertia coefficients in calculating the exciting fluid force 
through Morison's equation / 11_/, modified to take into 
account the motion of the structure. Their calculations have 
shown that a resonant motion is possible and that the dynamic 
deflection of a flexible structure can considerably exceed its 
static deflection. 

Verly and Every / 12_/ measured wave-induced stress on 

similar rigid and flexible vertical cylinders in a wave 

channel at relatively low Keulegan-Carpenter and Reynolds 

numbers. Even though they were unable to correlate their 

data with any suitable parameter governing the motion, they 

concluded that the vibration is caused by the cylinder's 

response to eddy shedding and that there is no fluid-structure 

interaction. They found that the vibration occurs if 0 /f D 

m n 
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(U represents the maximum velocity of flow or wave in a 
cycle; and D, the diameter of the cylinder) is greater than 
about unity for any natural frequency, wave frequency, and 
damping. The U m /f n D parameter never reached high enough 
values in Verly and Every' s experiments for the cylinder to 
undergo synchronized oscillations as the present investigation 
has shown . 

McConnell and Wilson / 13_/ conducted experiments similar 

to those of Laird / 8_/ by sinusoidally driving a circular 

cylinder in still water. These experiments were conducted 

at rather small Reynolds numbers (Re = U D/v) with 0.25 in and 

m 

0.50 in cylinders. Furthermore, these investigators measured 
only the amplitude of oscillation and calculated the lift 
coefficient at resonance from a slightly-modified version of 
an equation developed by Sarpkaya / 14_/. They have concluded 
that U m /f n D is one of the most important parameters governing 
the resonant oscillations. 

Sawaragi et al. / 15 / investigated the in-line and trans- 
verse dynamic response of a cantilevered circular cylinder, 
with a concentrated mass at its top, in waves of small ampli- 
tude. The Reynolds number ranged from 1,500 to 6,200 and the 
rms value of the Keulegan-Carpenter number (calculated over 
the submerged length of the cylinder through the use of rms 
value of the maximum of the horizontal velocities) ranged 
from 2 to about 20. They have approximated the lift coeffi- 
cient for a rigid cylinder by a Rayleigh distribution and 
calculated the dynamic response of the test cylinder. The 
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results are of limited value since the interaction between 
the synchronization and the force amplification was ignored 
and since no more than three vortices were shed during a half 
cycle. Furthermore, the rapid variation of the Reynolds 
number, Keulegan-Carpenter number, vortex-shedding coherence, 
and lift force with depth, in a range of Keulegan-Carpenter 
numbers where the phenomenon is least understood, obscured the 
role played by the governing parameters . 

It is clear from the foregoing that there are no well- 
established analytical methods either to analyze the data or 
to predict the characteristics of the flow or of the vibra- 
tion. Thus, it was decided to undertake a comprehensive in- 
vestigation of the dynamic response of elastically-mounted 
smooth and rough cylinders in harmonic flow particularly at 
high Reynolds numbers and Keulegan-Carpenter numbers. 
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Ill . EXPERIMENTAL EQUIPMENT AND PROCEDURES 



A. U-SHAPED OSCILLATING-FLOW TUNNEL 

The equipment used to generate the harmonically oscilla- 
ting flow has been extensively used at this facility over the 
past six years / 2_/, / 3_/. Only the salient features, most 
recent modifications, as well as the adaptation for this work, 
are briefly described in the following. 

The length of the U-shaped vertical water tunnel has been 
increased from 30 ft to 35 ft and its height from 16 ft to 24 
ft. Previously, a butterfly valve arrangement at the top of 
one of the legs was used to initiate the oscillations. During 
the past year, the tunnel has been modified so that oscilla- 
tions can be generated and maintained indefinitely at the 
desired amplitude. For this purpose the output of a 2 Hp fan 
was connected to the top of one of the legs of the tunnel with 
a large pipe. A small butterfly valve, placed in a special 
housing between the top of the tunnel and the supply line, 
oscillated harmonically at a frequency equal to the natural 
frequency of the oscillations in the tunnel. The oscillation 
of the valve was perfectly synchronized with that of the flow 
through the use of a feedback control system. The output of 
a pressure transducer (sensing the instantaneous acceleration 
of the flow) was connected to an electronic speed-control unit 
coupled to a DC motor oscillating the valve plate within 0.001 
seconds. The amplitude of oscillations was varied by 
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constructing or enlarging an orifice at the exit of the fan. 

The flow oscillated at a given amplitude as long as desired. 

Smooth and sand-roughened aluminum cylinders of 4 in 
diameter were used in the experiments . Each cylinder was 
rigidly fixed to two load cells and the low-friction linear 
bearings mounted on two matched helical springs. 

B. CYLINDER SUPPORTS AND HOUSINGS 

Initially, beam springs were considered for the elastic 
cylinder supports. However, the requirements for rigidity, 
strength and torsional stability in the in-line direction 
quickly led to the use of helical steel springs . In the formal 
design, the cylinder was rigidly fixed to two low-friction 
linear bearings resting on two matched helical springs. The 
motion of the bearings was constrained by two rods perpendicu- 
lar to the direction of the flow as shown in Fig. 1. Detailed 
drawings of the system are shown in Figs. 2, 3, and 4. 

The force transducers had a capacity of 500 lbs with an 
overload capacity of 200 percent. The transducers were mounted 
as shown in Fig. 4 to the aluminum split-bearing housing 
blocks. Each block was manufactured to accommodate an adjust- 
able ball bushing with a nominal 2.5 in bore and a length of 
5 in. The bearings rode on 2.5 in diameter stainless steel, 
hardened and ground, solid shafts. Adjustments of the bearing 
diameter eliminated bearing chatter. The springs were centered 
around the shafts by a retaining ring on the housing and by 
the yoke support at the bottom. 
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The yoke supports allowed the blocks and springs to be 
raised or lowered to center the cylinders vertically in the 
test section and permitted the use of various cylinders, 
springs and weights while accommodating the effects of the 
buoyant force. Each yoke was positioned by brass rods sus- 
pended through the top of the housing by another yoke which 
was positioned with precision by a threaded third brass rod 
passing through the upper yoke and lifting against the housing 
top. Alignment of mounts and rods was accomplished by a jig 
which was placed between the two blocks, leveled, then raised, 
and lowered to position the two stainless steel rods. 

Figure 5 shows a dimensioned drawing of one of the two 
support housings. Each was welded of 3/8 in aluminum plate. 
The windows were covered by removable plexiglass plates for 
access for calibration and maintenance, or, when mounted in 
place, for the observation of cylinder movement during tunnel 
operation. Vents and drains were installed at top and bottom 
of the housing to remove air or water in filling or draining. 
The housings were flanged and bolted to the U- tunnel wall with 
O-ring seals for water proofing as may be seen in Fig. 3. 

Linear-variable-displacement transducers (LVDT) were in- 
stalled between each bearing block and housing. Each LVDT 
was energized by a six volt DC regulated power supply. The 
LVDT's provided displacement from the equilibrium position and 
also helped to verify the simultaneous and in-phase vertical 
oscillation of both ends of the cylinder. 
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The circular cylinder models used for this study were turned 
on a lathe from aluminum pipes and polished to a mirror shine. 
The length of each cylinder was such as to allow 1/32 in 
clearance between the tunnel wall and each end. As part of the 
total length, smooth circular flanges 6.0 in diameter and 1/8 in 
thick were attached to the ends of the cylinder to eliminate 
the end effects. 

All cylinders were hollow, and each contained a beam type 
accelerometer. Figure 6 shows a cutaway drawing of a cylinder 
with the accelerometer installed in the exact center of the 
cylinder. The electrical connection to the accelerometer is 
led through a notch in the male fitting at the end of the 
cylinder. Sealing of the unit is accomplished by O-rings and 
threaded fittings. 

A 6.0 in circular access hole was cut in the tunnel walls 
at the test section. Each hole was covered with a plate of 
identical size. A slot of 1.0 in by 5.0 in was cut in each 
cover. The covers allowed the installation of the cylinders, 
and the slots allowed the unobstructed linear oscillation of 
the force-transducer arms attached to the cylinder. 

C. MEASUREMENT OF BASIC PARAMETERS 

For the determination of k, a known weight was suspended 
from the center of the cylinder, and the displacement was 
measured. The natural frequency in air, f , and the natural 
frequency in water, f , were determined by plucking the 
cylinder and allowing it to oscillate freely in the particular 
fluid medium. The displacement was recorded on one channel 
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recorder. A sample trace is shown in Fig. 7. Along with the 
frequencies , the logarithmic decrement and the damping factors , 

5 , and, were determined in air and water, respectively. 

The mass, M, of the system was then calculated from 

M = k /(2irf ) 2 (1) 

lid 

and the added mass of the system was determined from 

M + M = k/ ( 2 7T f ) 2 = M(f /f ) 2 (2) 

a nw na nw 

For the purpose of verification of the measured values, 
for the cylinder was computed through the use of the observed 
values of frequencies and found to be within one tenth of one 
percent of that given by the potential flow theory and inde- 
pendent experiments . * 

A trace was also made of the output of the accelerometer 
oscillating in free vibration. The frequency of this device 
was found to be 44.5 Hz with a damping factor of ? = 0.0022. 
According to Thompson /~16 /, for frequencies up to 5 Hz this 
accelerometer would provide true force readings with a magni- 
fication factor nearly equal to unity. The frequency of oscil- 
lation in the present investigation ranged from 1.5 Hz to 4.0 Hz. 



* The theoretical value of the added mass coefficient for 
a circular cylinder is equal to 1 . 0 . 
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Amplitude, acceleration, elevation or velocity of the fluid 
was a matter of interpretation of the signal received from the 
same instrument. In this experiment use was made of a differ- 
ential pressure transducer with pressure taps located on the 
two legs of the tunnel. Knowing the height of the mean water 
level above the taps, H, the above values were solved for by 
application of Bernoulli's equation. For example, the virtual 
amplitude is given by / 2 / 



A 



m 




AP 

1 - — ( 2tt/T) 2 H 

g 



(3) 



in which the gravitational acceleration, g, and the period of 
flow oscillation, T, remain constant. 

By examination of the free body diagram shown in Fig. 8, 
it is seen that the only transverse forces acting on the 
cylinder are those reaction forces from the two transducers, 
the fluid force and the product of the cylinder mass and 
acceleration. Therefore, the fluid force is given by 

F^ = F^ + - My (4) 

The sensitivity of the accelerometer was adjusted to yield 
F = 0 while the cylinder freely oscillated in air. This 

J-J 
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method allowed the measurement of the net fluid force with an 
error less than 1.5 percent. 

D. ELECTRONIC CIRCUITRY 

The output of each of the two LVDT's, after amplification 
by a low-gain amplifier, was displayed on an eight-channel 
recorder as in Fig. 9. One of the two signals was additionally 
displayed on a three-channel recorder. 

The output of the force transducers was separately ampli- 
fied by a carrier amplifier, filtered by a Krohn-Hite elec- 
tronic filter and displayed on the eight-channel recorder. 

Then the two signals were added together by an electronic 
summing circuit, amplified again by a low level preamplifier 
and displayed on the eight-channel recorder. 

The accelerometer output was similarly amplified by a 
carrier amplifier and filtered by a low pass filter. This 
signal was also displayed on the eight-channel recorder. It 
was then subtracted from the Siam of the force signals through 
the use of a summing-differencing electronic circuit. This 
resultant signal, representing the fluid force on the cylin- 
der, , was again amplified by a low level preamplifier and 
displayed on both the three- and eight-channel recorders. 

Figure 10 is a block diagram of the circuits described above. 

The signal from the differential pressure transducer was 
of such quality that it did not require any filtering. As 
shown in Fig. 11, it was amplified by a carrier amplifier and 
displayed on the two aforementioned recorders. 
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E. PROCEDURES 



Procedures for a series of runs were as follows: 

1. The cylinder was aligned with the tunnel by adjusting 
the bearings. The sensitivities of the force transducers 
were matched by hanging a known weight from the exact center 
of the cylinder. At the same time a calibration factor was 
determined by measuring the displacement and comparing it 
with the output of the LVDT of the three-channel recorder. 

2. The accelerometer was calibrated by plucking the 
cylinder (quickly releasing the cylinder from a displaced 
position and allowing it to oscillate freely) and adjusting 
the amplifier gain until the fluid force in the air was 
nulled to within 1. 5 percent deviation. 

3. f and £ were also determined by plucking the 

na a 1 c 3 

cylinder in air . 

4. The weight was then suspended from the cylinder center, 
this time for the measurement of a calibration factor in 
lbf/mm of the three-channel-recorder output at a particular 
gain setting. 

5. The U- tunnel and cylinder-support housings were again 
filled to the top of the test section, and all air was vented 
from the housings. 

6. The cylinder was then lowered to the center of the 
test section by simultaneously lowering both sets of spring 
yokes, each the same number of turns, to compensate for the 
buoyant force on the displaced cylinder. 

7. The cylinder was plucked, and f nw and were calcu- 
lated from the free-oscillation traces. 
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8. The access to the horizontal test section above the 



cylinder, through which the weights were lowered to the 
cylinder, was closed. The tunnel was filled with water from 
a sump tank located near the apparatus . Filling was acceler- 
ated by a large submerged sump pump. 

9. The differential-pressure transducer was bled of any 
entrained air and its calibration verified. 

10. The tunnel oscillations were started by actuating the 
control system and the blower. 

11. The amplitude of flow oscillations were set at a 
desired value and the forces, displacements, and flow acceler- 
ation were recorded simultaneously. 

12. Three or more runs were made at a paper speed of 25 
ram/sec at each combination of spring, cylinder and flow 
amplitude . 

13 . At the end of each set of runs , f and c, were 

nw w 

measured to verify the initial values before draining the 
tunnel . 

Two additional runs were made with the cylinder mounts' 
fixed to verify the results previously obtained by Sarpkaya 
/ 2 J , The force transducers were turned 90 degrees and 
pinned to measure the in-line forces. Also readings of these 
forces were taken with the cylinder free to oscillate. Data 
from these runs proved to be as expected, and they were not 
repeated for the larger cylinder. 
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F. DATA REDUCTION 



Test runs were selected for evaluation on the basis of 

oscillation amplitude, Reynolds number, and Keulegan-Carpenter 

number. For each run evaluated, the peak lift force, F T „ ; 

LM 

peak displacement, Y^ ; shortest force- and displacement- 
periods, T v and T c ; and flow amplitude were recorded for 
each half cycle of flow oscillation. Also recorded was such 
relevant information as calibration factors, cylinder dia- 
meters, gain settings, k, . f , etc. A computer program 
calculated the peak lift coefficient, C , as follows: 



c 



LM 



1 

2 



LM 



pU 2 ld 

m 



( 5 ) 



where p is the density of the fluid; and L , the length of 

the cylinder. Other variables computed included: A m , 

U , K , Y /D , U /f D , and the Strouhal number St. 
m M nr n 
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IV. ANALYSIS 



The amplitudes of the lift coefficient and the relative 
displacement may be shown to be governed by 

tC LM' V D> = £ i {0 m T/D ' d2/vT ' V £ m D ' V D ' V M / pL ° 2} (6) 

where 



C T = F_ M /(0.5pLDU^) (7) 

LM LM m 

and represents the amplitude of velocity oscillations; 
f , the in-water natural frequency of the oscillating system; 
D, the diameter of the cylinder; v , the kinematic viscosity; 
T, the period of the flow oscillations; k g , the roughness 
height, C , the damping ratio in air; M, the mass of the 

3 

oscillating system; p , the density of the fluid; and L, the 

length of the cylinder. The last two parameters in Eq. (6) 

2 

may be combined as M? /(pLD ) on the basis of experience with 

3 

hydroelastic oscillations in steady flow / — 7 

The lift coefficient and the Strouhal number for a non- 
oscillating cylinder are governed by 

{S ° ' C LM } = g i {U m T/D ' d2 / vT ' k s /D} (8) 
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In fact, the three parameters in Eq. (8) determine uniquely all 
the characteristics of the base flow. Missing from the above 
relationship is the effect of the correlation length. For 
smooth circular cylinders the lock-in or the forced-oscillatory 
motion is accompanied by increased coherence in the shedding 
along the body. This increased correlation was thought to be 
primarily responsible in increasing the force and hence the 

amplitude to some limit set by the flow. Experiments with 
bodies with fixed seperation lines (and suitable after body) 
have shown / 17 / that the increase in correlation length can 
account for only a small fraction of the large increase in the 
transverse force near lock-in. The fluid response to the body 
motion or the interaction between the fluid motion and the 
body motion is the primary reason for the lock-in and the in- 
crease in the transverse force. 

The problem under consideration requires the solution of 
the separated, time dependent flow about an oscillating cylin- 
der. Such a solution is hard to come by in view of the fact 
that a solution does not exist even for steady flow about a 
cylinder at rest. Thus, a phenomenologically based analysis 
coupled with experimental observations is the best one can 
perform at present. Hopefully, the insight gained through 
such an analysis will eventually lead to the formulation of 
an heuristic model. 

The equation of motion for an elastically-mounted rigid 
cylinder of mass M, damping coefficient £, circular natural 
frequency w , and spring constant k, is given by 
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My + 2M£<i^y + ky = F(t) 



(9) 



where y represents the transverse displacement of the cylinder 

and F ( t) the exciting force or the lift force. In general 

F ( t) depends on the parameters given by Eq. (6) or simply 

2 

U r = U m /f n D and M?/(pLD ). For sake of simplicity, F(t) may be 
assumed to be given by 

F ( t) = 0.5pLDU 2 C TM sin (co t + 9) (10) 

m LM v 

in which ai v = 2 irf and 9 is the phase angle. However, even 
during the period of nearly perfect synchronization, the lift 
force does not retain a constant amplitude and frequency. 
Furthermore, the velocity U varies as U = U m sin 2irt/T during 
a given cycle. In general, it is necessary to consider the 
harmonics of the lift force together with their spectral 
energy. Prior to such a detailed analysis, it will be hy- 
pothesized and subsequently verified by the present results 
that the rate of energy transfer from the fluid to the oscil- 
lating body is proportional to the rms value of the velocity. 
Writing , 

2 

T T U 

i f U 2 dt = | / U 2 sin 2 2irf t dt = -f (11) 

i 0 1 o 

and combining with Eq. (10) , one has 

F ( t) = j pLDU 2 C LM sin(w v t + 9) (12) 
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The equation of motion may be written as 

My + cy + ky = F_,.sin(to t + 9) 

LM v 

The steady state solution of this equation is given by 

y = A sin (to t + 9 ) + B cos (to t + 0) 

Substituting Eq . (14) in Eq. (13) and solving for A and 

one has 



A = F 



(k - Mto ) 
v 



LM ,, ,,2,2 , 2 2 

(k - Mto ) + to c 

v v 



and 



B = F 



to c 

V 



LM Tv „ 2, 2 , 22 

(k - Mto ) + to c 

v v 



Then, Eq. (14) may be written as 



y = Y M cos{(to v t + 0) - ip) } 



in which 



M 



- / 



2 2 
A + B 



( 13 ) 



(14) 

B 



(15) 



(16) 



(17) 



(18) 
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and 



♦ = 



tan ^ 



0) 



25 



v 



fcO 



n 



03 



1 - (- 



03 



^> 2 



n 



(19) 



and 



y m - 



F LM //k 



( 20 ) 



/ 



0) . 



03 



a - ( -^-> 2 > 2 + <2c” v 



n 



03 



n 



then Eq. (17) reduces to 



Y = 



P LM^ 



cos (w t + 0 

V 



- 1(0 



/ 



03 n « 03 0 

(1 - ( " 77 ~ — ) 2 } 2 + ( 2 Z -^-) 2 



03 



n 



03 



n 



Noting that k = M(27Tf n ) , Eq. (20) reduces to 



- y m 
y m - IT 



( V f n D) 2 C LM pL ° 2 



,2„ /" r , , “v ,2,2 . “v ,2 



( 22 ) 



16tt M / {1 - ( — — ) ) + (2? 



03 



n 



n 



At resonance 03 ~o3 , and Eq. (22) becomes 
v n ^ 



- y m 
y m d 



(U /f D) C TM 
m n LM 



32tt 2 M?/pLD 2 



(23) 



( 21 ) 
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Equation (22) may also be written as 



M 



'LM 



i /- 2., / hi _ » CO ~ 

16 ” M / U - <-V> 2 + (2t^) 2 

If n n 



pLD (O n /f n D) 



(24) 



The above equation is used to compare the measured and calcu- 
lated values of the lift coefficient in terms of the known 
and measured parameters appearing on the right-hand side of 
the equation. 

The foregoing comparisons between the measured and calcu- 
lated values of Y../D and C T „. are only for the peak values of 

M LM - 1 

the said parameters. Evidently, it is equally important to 
compare the evolution of the displacement as a function of 
time using the measured force as the input function. 

Among the various techniques commonly used in the vibration 
analysis, Duhamel's superposition principle enables one to 
predict the displacement of a linearly-damped and elastically- 
mounted cylinder through the use of the following integral 
/~16 7- 



y = 



(y — i 

v o Mui 



F( + )e 



-5Vt-V 



n o 



r 



sin / 1 - \ 



co t dt ) cos to t 
n v v n 



1-5 



n no 



-5u <+-+ ) 
T F ( t ) e n v 



/ 



cosJ 1- c co t d+ )sinoo f 
n v v n 



1 - 



Equation (25) is used to compare the predicted and measured 
relative displacements for a series of selected runs. 



( 25 ) 
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V. DISCUSSION OF RESULTS 



The results will be discussed first in terms of the maxi- 
mum values of the lift coefficient and the amplitude of oscil- 
lation. 

Figures 12 through 14 show the variation of with re- 

spect to the Keulegan-Carpenter number. The variation of the 
same parameter with respect to the reduced velocity U r is 
shown in Figures 15 through 17 for a smooth cylinder and for 
a rough cylinder with two different damping parameters. All 
of these figures show that decreases rapidly with in- 
creasing K and U^. . Furthermore, is considerably larger 

for the rough cylinder. The change in damping does not change 
the value of as evidenced by a comparison of Figures 13 

and 14 and of Figures 16 and 17. The smooth cylinder has not 
been tested with different damping ratios. However, it is 
expected that the data shown in Figures 12 and 15 will not 
vary with damping . 

The variation of with the ratio of the vortex shedding 

frequency to the cylinder oscillation frequency is shown in 
Figures 18 and 19 for a smooth and a rough cylinder. Evidently, 
as vortex shedding frequency becomes nearly equal the cylinder 
oscillation frequency, the lift coefficient increases drama- 
tically. It is this particular nature of the vortex- induced 
oscillations that has been referred to in the literature as 
self-excitation. A vortex-excited oscillation is actually a 
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forced one having a self-excited character also to some degree 
due to lift-force amplification through nonlinear interactions. 
By definition, self-excited oscillation is one where the 
alternating force that sustains the motion is created or con- 
trolled by the motion itself; when the motion stops, the al- 
ternating force disappears. For the phenomenon discussed 
herein the alternating lift force does not disappear when there 
is no oscillation. In fact, when there is no alternating force 
(a flat plate normal to the flow) there is no lock-in. Appar- 
ently, the often-used definition of "self-excited oscillation" 
is a misnomer. 

The calculated values of in the synchronization region, 

through the use of Eq. (24) and the measured values of Y^/D 
and £ , are compared with those measured directly in Figure 20. 
Clearly, the measured and predicted values are in good agree- 
ments in the region where Eq. (24) is valid. In reality, a 
number of frequencies exists for the vortex shedding even in 
the synchronization region (0.95 < f v /f < 1.05) . This is 
particularly true for smooth cylinders where the vortex 
shedding along the cylinder is not well correlated. For rough 
cylinders the correlation along the length of the cylinder is 
almost perfect / 2_/. This is partly the reason for the 
better agreement between the measured and calculated values 

of C T „ shown in Figure 20. 

LM 

The relative displacement of the cylinder is examined 
for both smooth and rough cylinders in terms of the reduced 
velocity, Keulegan-Carpenter number and the Reynolds number. 
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Figures 21 through 23 show the variation of Y^/D as a function 
of U r . The maximum amplitude at synchronization is reached 
at about U =5. 9 for the smooth cylinder and at about U r =5.6 
for the rough cylinders. Experiments with steady flow about a 
tranversely-oscillating smooth cylinder have shown that / 1 _/ 
perfect synchronization takes place for values between 5 
and 7. Evidently, synchronization in harmonically oscillating 
flows falls in the lower range of U r values just cited. 

The variation of Y M /D with the Keulegan-Carpenter number 
is shown in Figures 24 through 26. Apparently, the synchroni- 
zation for the 4 in cylinder occurs at a Keulegan-Carpenter 
number of about 50 . This result is not unexpected as will be 
demonstrated below. 

The Strouhal number for a harmonically oscillating flow 
may be written as 

f D 

St = -g— (28) 

m 

multiplying and dividing the right-hand side of this equation 
with T (the period of oscillations of the water in the tunnel) 
one has 



St 



f v T 
K 



(29) 



at synchronization f is nearly equal to f n » For the 4 in 
cylinder f n = 1.6 and T = 6. sec. With an average Strouhal 
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/ 



number of 0.2 Eq. (29) yields K = 48 which compares quite well 
with the K value obtained from Figures 24-26. 

The variation of Y^/D with the Reynolds number is shown 
in Figures 27 through 29. For the particular cylinder the 
synchronization occurs at a Reynolds number about 100,000. A 
comparison of Figures 24-26 with 27-29 shows that the depen- 
dence of y m /D on K is exactly the same as that of Y^/D on the 
Reynolds number. Even though, experiments have not been per- 
formed at a given K for different Reynolds numbers the results 
with steady flow experiments show that / 7 _/ Y M /D at the 
synchronization does not depend on the Reynolds number in the 
range of Reynolds numbers from about 1,000 to 10 ^ / 7_/. Thus, 
Figures 27-29 actually reflect the variation of Y^/D with the 
Keulegan-Carpenter number. 

The amplification factor, defined by Y^k/F , in the syn- 
chronization region is plotted in Figures 30-32 as a function 
of f^f for both smooth and rough cylinders. Considering 
the difficulty of establishing and maintaining synchronization, 
the experimental results are in excellent agreement with that 
predicted by Eq. (20) . 

As noted earlier, synchronization is possible only when 
the frequency of the exciting lift force is close to the 
natural frequency of the oscillating system. Consequently, 
the Strouhal number plays a major role in the phenomenon. 
Figures 33 through 35 show the variations of Y M /D and St with 
U . For the smooth cylinder the Strouhal number decreases 
from 0.2 to about 0.16 as one approaches synchronization. 
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Thereafter, the Strouhal number gradually increases to about 
0.2. The variation of the vortex shedding frequency and 
hence the Strouhal number for smooth cylinders is strongly 
effected by the change of the correlation length / 7_/. On 
the other hand, roughness provides perfect correlation and 
Strouhal number remains nearly constant at St 1 0.2 as seen 
in Figures 34 and 35. Apparently, roughness enables one to 
eliminate from consideration the effect of the correlation 
length. Initially experiments were conducted with a smooth 
cylinder. It was supposed that roughness will complicate the 
matters further. Contrary to this supposition, roughness and 
the harmonic nature of the flow further delineated the most 
important parameters governing the hydroelastic oscillations 
of cylinders. 

Finally, Figure 36 shows the comparison of the instantan- 
eous values of the measured cylinder displacement with those 
predicted through the use of the Duhamel's integral / see 
Eq. (27)_/. The measured force and other physical charac- 
teristics of the oscillating system (damping, mass and natural 
frequency) were used as input in the Duhamel's integral. The 
close agreement between the measured and calculated dis- 
placements is not surprising. It simply proves that the 
various parameters such as mass, damping ratio, displacement, 
etc . , have been properly measured . 
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VI. CONCLUSIONS 



The theoretical and experimental investigation of the 
hydroelastic oscillations of a smooth and rough circular 
cylinder in harmonic flow warranted the following conclusions: 

1. An elastically-mounted cylinder may undergo synchronized 

oscillations when the reduced velocity U = U /f D is 

1 r m n 

equal to about 5.7. 

2. The synchronized oscillations for a smooth cylinder 
occur at an average Strouhal number of about 0.16. 

The Strouhal number for the rough cylinder remains 
constant at a value of about 0.2. 

3. In the region of synchronous oscillations, the lift 
coefficient is amplified by a factor of about 2 rela- 
tive to that for a fixed cylinder in the same flow. 

4. The analysis relating the governing parameters predicts 
fairly accurately the maximum amplitude of the oscil- 
lations for both the smooth and the rough cylinders in 
the synchronous region. 

5. The measured instantaneous values of the cylinder 
displacement are in excellent agreement with those 
predicted through the use of the Duhamel's integral. 
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TOP VIEW 





Fig. 1 Preliminary sketch of the supports and housings. 
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Fig. 2 



Side view of a cylinder-support mount 
and housings. (Dimensions in cm). 
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L. 
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Fig. 3 Top view of a cylinder-support mount and housing. 
(Dimensions in cm) . 
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Fig. 4 Side Cutaway view of a Load cell 
bearing housing. (Dimensions in 




Fig. 5 Support housing. 

(Dimensions in cm) . 
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Fig. 6 A cutaway sketch of a cylinder with 
accelerometer. 
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Fig. 7 Displacement trace of a "plucked" cylinder. 
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Fig. 8 Freebody diagram of the 
test cylinder. 
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Fig. 9 LVDT ( Linear-variable-displacement transducer) electronic circuit. 




Fig. 10 Force and accelerometer electronic 
circuits 
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Fig. 11 Differential pressure transducer 
and circuitry. 



50 



61 'II 



J L 

co 



a 

cm 



m 

XL 



(I 



CO 

o 

II 

< 



£00 

w 

e 
£0 
£ 0 © 
0 O 



a 

a 



o 

CO 



□ 



©0 



©0 



00 0 
0 



D0D 



a 

CD 



X 

z: 

ZD 

II 

DZ 



G @)0 
G 00 

0 © 

00 0 ® 

0 ©00 



0 00 
Q€ b 0 

Q® 0 9 0 



J L 



a 
; c\i 



CM 



Q3ynSU3N-W13 



51 



Fig. 12 C versus K for a smooth cylinder. 
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Fig. 13 C versus K for a rough cylinder with c =0.060 
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Fig. 14 C versus K for a rough cylinder with £ =0.03 
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Fig. 15 C T versus U for a smooth cylinder. 
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Fig. 16 C versus U for a rough cylinder with L =0.060 
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Fig. 17 C versus U for a rough cylinder with c, =0.03 
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Fig. 18 C versus _v for a smooth cylinder 
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Fig. 19 C versus f /f for a rough cylinder with t =0.060 
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Fig. 20 C T __ versus C T __ for smooth and 
17 LM LM 

rough cylinders in the synchronization 

region 



59 



1852 



CD 

t — 1 


i — 1 


■ 


CO 


t — 1 


• 


! 1 


o 


il 


II 


H, 


< 


ro 

C 

Oh 





L. 

co 

o 



J 1 I 

CD ZT 

O O 

0/MA 



ee 



OJ 



0 0 
0 0 
3)0 
€0 0 



_ O 



00 



0)0 

£0 



CO 



a 

x 

Ll_ 



€E0 

00 

0^0 

h-&) 

00) 

®SD 



CD 



0 



€D . 



OJ 



OJ 

o 



o 

o 



60 



Fig. 21 Y /D versus U for a smooth cylinder 
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Fig. 22 Y /D versus U for a rough cylinder with c =0.060 
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Fig. 23 Y /D versus U for a rough cylinder with C =0.030 
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Fig. 24 Y / D versus K for a smooth cylinder 
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Fig. 25 Y /D versus K for a rough cylinder with C =0.060 
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Fig. 26 Y /D versus K for a rough cylinder with C =0.03 
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Fig. 27 Y /D versus Re for a smooth cylinder 
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Fig. 28 Y /D versus Re for a rough cylinder with 1 =0.060 
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Fig. 29 Y /D versus Re for a rough cylinder with L =0.03 
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Fig. 30 Y..k/F T ,. versus f /f for a smooth 
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cylinder in the synchronization region 
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Fig. 31 Y^k/F^ versus f y /f n ^ or a rou< ?k cylinder 

with ^=0.060 in the synchronization region 



70 



YM*k/FLM 




FV/FN 



Fig. 32 versus f v /f n ^ or a roU( ?^ 

cylinder with C.=0.03 in the 
synchronization region 
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Fig. 33 Y /D and St versus U for a smooth cylinder 
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Fig. 34 Y /D and St versus U for a rough cylinder with C =0.060 
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Fig. 35 Y / D and St versus U for a rough cylinder with £ =0.03 
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